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We show that discotics, lying deep in the columnar phase, can exhibit an x-ray scattering pattern 
which mimics that of a somewhat unusual smectic liquid crystal. This exotic, new glassy phase 
of columnar liquid crystals, which we call a "hybrid columnar Bragg glass", can be achieved by 
confining a columnar liquid crystal in an anisotropic random environment of e.g., strained aerogel. 
Long-ranged orientational order in this phase makes single domain x-ray scattering possible, from 
which a wealth of information could be extracted. We give detailed quantitative predictions for the 
scattering pattern in addition to exponents characterizing anomalous elasticity of the system. 
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Until now, the x-ray scattering pattern given in Fig.^ 
would be identified with a system in a somewhat unusual 
smectic phase with short-ranged translational order and 
long-ranged orientational order within the smectic layers; 
i.e., a smectic composed of nematic, rather than liquid, 
layers. The set of on-axis quasi-sharp Bragg peaks along 
qh is a signature of the quasi-long-ranged translational or- 
der (i.e., the periodicity of the layering) that is character- 
istic of the bulk smectic phase. The presence of the other, 
broadened, peaks and the azimuthal anisotropy about the 
qh axis respectively indicate the incipient short-ranged 
translational order and the long-ranged orientational or- 
der within the smectic layers oriented perpendicular to 

qh - q. 
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FIG. 1. X-ray scattering pattern in the L plane for a class 
of hybrid columnar Bragg glass. 

In this Letter we predict the existence of a remarkable 
new "hybrid columnar Bragg glass" phase (HCBG) 
which, despite differing fundamentally from the smectic 
phase described above, shares the same qualitative scat- 
tering pattern illustrated in Fig|l]. Such mimickry of one 
phase by a completely different other phase is unprece- 
dented. 

Columnar phases in pure, bulk (i.e., quenched-disorder 
free) liquid crystals are phases that have long-ranged 
translational order in two-directions, and short-ranged 
translational order (i.e., liquid-like correlations) in the 
third. I.e., they are regular two-dimensional lattices of 



one-dimensional liquid columns (Fig.^|(a)). In this Let- 
ter, we show that when such a system is confined in an 
anisotropic quenched random environment, e.g., strained 
aerogel ||, it becomes translationally disordered, but re- 
mains topologically ordered (i.e., free of topological de- 
fects like dislocations). This novel state is the HCBG. 

Like the smectic phase, the HCBG has translational or- 
der that is quasi-long-ranged in one direction and short- 
ranged in another, as implied by Fig.|l|. As illustrated 
in Fig.||(b), the columns remain in roughly equidistant 
rows perpendicular to the stretching direction, but lose 
long-ranged translational order within each row. Nev- 
ertheless, the hexagonal orientational order is preserved, 
albeit uniaxially distorted due to the strectch. However, 
there are a number of fundamental differences between 
the two phases. First, unlike the smectic phase, the ex- 
ponent rye characterizing the shape of the quasi-long- 
ranged translational order peaks at G, which are given 
by /(q) cx |q — G|~ 3+,,G , is independent of tempera- 
ture. This property would allow the two phases' scat- 
tering patterns to be distinguished through comparison 
of the lineshapes as temperature is varied. Secondly, the 
correlations of the quasi-long-ranged order in the HCBG 
scale isotropically in space, in contrast to the well-known 
strongly anisotropic scaling of these correlations in the 
smectic. The third and most crucial difference between 
the two phases is their topolological order, which distin- 
guishes their elasticities but not their scattering. Specif- 
ically, the absence of the extra direction of long-ranged 
translational order, in the smectic phase, is caused by free 
topological defects, namely unbound dislocations with 
Burgers vectors along the smectic planes. Although the 
HCBG also exhibits translational order that is quasi- 
long-ranged in one direction and short-ranged in another, 
at long length scales it is distinguished from the smectic 
by being free of these unbound dislocations, which would 
otherwise destroy the columnar phase topology of HCBG. 
One important experimental consequence of this absence 
of free dislocations is that the HCBG retains elastic re- 
sistance to distortions in the extra direction, albeit, as 
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we discuss below, of a very strange, anomalous sort. 

Of course, for sufficiently strong disorder free disloca- 
tions will eventually proliferate. The anisotropy, imposed 
by the strained aerogel, leads to the interesting possibil- 
ity that dislocations with their Burgers vectors in the soft 
direction may unbind (with increasing disorder) before 
those with their Burgers vectors in the hard direction, 
leading to the sequence of disorder-driven phase transi- 
tions : HCBG — > m=l smectic Bragg glass B — * nematic 
elastic glass || with increasing aerogel density. 

The rest of this Letter gives a more detailed theoretical 
description of the HCBG phase, including x-ray correla- 
tion lengths and universal exponents characterizing the 
anomalous elasticity. We relegate the technical details to 
a future publication || . 

Our model for a columnar phase consists of disc-shaped 
molecules with normals aligned along the z direction. 
The discs form a hexagonal lattice in the xy (_L) plane 
and have liquid like correlations along z as illustrated in 
Fig.^(a). We assume, and verify a posteriori, that despite 
considerable distortion, for sufficiently weak quenched 
disorder, the columnar phase topology is stable, i.e. our 
discotic liquid crystal remains free of unbound disloca- 
tions. Consequently this system can be described within 
an elastic theory, with a two-component (x and y) lat- 
tice site displacement vector u(r) and the discotic direc- 
tor n(r) (the normal to the discs) as the only important 
long length-scale degrees of freedom. The disordering 
tendency of the aerogel is two-fold: the strands act both 
to randomly pin the columnar lattice (u(r)) and to dis- 
tort the orientations of the disc normals (n(r)). Our 
starting Hamiltonian is that of a pure hexagonal discotic 
in isotropic aerogel 
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\d z u - 5n\ 2 + -Xu^ufj + im±u± 



^e^^(r)e* G - u «-(g(r).n) 



H F [h}. (1) 



where uj^ — \{d^Uj + dfui — dkUid^Uj) is the rotation- 
ally invariant symmetric strain tensor, <5n(r) = n(r) — z, 
the B± term reflects the tendency of the molecular di- 
rector (disc normal) n(r) to lie along the local tangent 
t w z + d z u to the liquid-like columns, Hp[h) is the 
Frank free energy of the molecular directors, and Vi(r) 
is a complex random pinning potential that couples to 
lattice site fluctuations along the reciprocal lattice basis 
vector, Gi. At long length-scales its correlations can be 
accurately represented as zero-mean with short-ranged, 
Gaussian statistics: Vi(r)V*{r') = A v 6±S d (r ~ r'), |j 
throughout this paper x denotes a quenched average over 
the disorder of the quantity x, while (x) denotes a ther- 
mal average. The last term describes the tendency of 
the disc normals n(r) to align along the random local 
aerogel strand directed along g(r). This "random tilt" 
disorder is described by short-ranged, isotropic correla- 



tions 5i (r) 3j (r') = l/2\/AM d (r - ""') @. Ay and A 
are phenomenological parameters which, in the simplest 
microscopic model, are proportional to the aerogel den- 
sity only, pa- As for smectics H, these two types of 
disorder have important long distance effects. 




Sir'. 



r 



% 



(a) 



- stretch 
(b) 



FIG. 2. (a) Alignment of the lattice for a uniaxial stretch, 
(b) Schematic of the distorted lattice in real space. 

A detailed analysis |||| has shown that this system 
exhibits a columnar Bragg glass phase with only short- 
ranged translational order. However, this changes if the 
aerogel is anisotropic. Aerogel anisotropy could be real- 
ized e.g., by applying a strain to the strands. For het- 
erotropic alignment between the disc normals and strands 
(assumed throughout) a uniaxial compression will lead 
to a phase in the same universality class as the Bragg 
glass phase of an Abrikosov flux lattice |I]J^] , with quasi- 
long-ranged translational order in both directions of the 
_L plane §. The more interesting HCBG phase with 
quasi-long-ranged order in only one _L direction can be 
obtained by applying a uniaxial stretch to the strands. 
For homeotropic alignment of strands and disc normals 
the two phases reverse with respect to stretch and com- 
pression leaving all of our other predictions unchanged. 

Uniaxial stretch of the aerogel strands (along eh causes 
the disc normals to align _L to the axis of stretch (FigJ|). 
This can be accounted for by the addition of a term 
Jr 7(eh-n) 2 to the Hamiltonian, where 7 is a phenomeno- 
logical parameter, which we assume to be a monotoni- 
cally increasing function of pa and strain. The uniax- 
ial stretch also breaks both the hexagonal symmetry of 
the lattice, and its rotation invariance. Thus, we are 
forced to consider both a more general lattice structure 
and an elastic Hamiltonian that need not, and indeed, 
will not be invariant under rotations of the lattice. We 
can take into account both of these effects by replac- 
ing the harmonic pieces of the hexagonal elastic energy 
§ r {\^ u ii u fj + ^ u ij u ij) with the more general harmonic 
elastic energy | J Cijkid^Ujd^ui, where the elastic con- 
stant tensor Cijki is not symmetric under interchange of 
its first two or second two indices, due to the lack of in- 
plane rotation invariance just discussed. The terms cubic 
and quartic in u in Eq.[l] must also be so generalized of 
course. Because of the B± term in Eq.[l], fluctuations of n 
from the local column tangent t are small, i.e., i5n « d z u. 
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We have analysed |6) this generalization of the model, 
Eq.|l], using renormalization group (RG) methods ||. 
One of the most surprising conclusions of this analysis is 
that, at long length scales, fluctuations Uh = ■ u along 
the direction of stretch decouple from those u s = e s • u 
orthogonal to this direction, where we have denoted the 
axis of stretch, e^, as "hard"(h) and called the other _L 
axis, orthogonal to ll soft"(s), i.e. = (r^,r s ). That 
is, all couplings between u s and that are present in 
the full elastic tensor Cijki flow to zero upon renormaliza- 
tion, leaving C ljk i in the form C ljk i = B ss S is Sj s 5k s Si s + 

BshSisSjhSksfilh + BhsSihSjsSkhSls + BhhSihSjhSkhSlh- 

The total Hamiltonian for the system can therefore be 
expressed as a sum of decoupled Hamiltonians for Uh and 
u s : H tot [u h ,u s ] = H X Y[uh] + H m=1 [u s ], with 



H XY = - 



l\d z u h \ 2 + B sh \d s u h \ 2 + B hh \d h u h \ 



(Gi-e h )u h (r) 



H m =i = J y(d 2 u s ) 2 + B ss (d s u s - ^(d z u s )'< 



+ B hs (d h u s ) 2 - g z (g- e s )(d z u s ) 



(2a) 



(2b) 



These two parts of H to t, which describe the fluctuations 
of the "hard" and "soft" phonon fields Uh and u s , are 
not new. Hamiltonians of precisely this form have been 
previously used to describe the random field XY-model 
g and the "m=l smectic Bragg glass" j|, respectively, 
and have been studied extensively. However, a colum- 
nar phase confined in anisotropic aerogel, whose Hamil- 
tonian, H to t, is a combination, or hybrid °f the two, is 
entirely novel. 

The Uh fluctuations of our system are the same 
as those of a slightly anisotropic (although isotropic 
in scaling) random field XY-model and are given by 
(K(r) - U h (0)) 2 ) = C(d)lnr/G 2 Qh , where G oh is the 
lattice spacing of the projection of the discotic recip- 
rocal lattice onto the hard axis. They diverge loga- 
rithmically as a function of distance, implying that the 
translational order along the hard direction is quasi-long- 
ranged. While these elastic distortions are reminiscent of 
the famous Landau-Peierls In r fluctuations of bulk smec- 
tics they differ crucially in two ways. Firstly, they are 
disorder, rather than thermally, driven with C(d) uni- 
versal (C(3) pa f.f) and the logarithm persisting in all 
2 < d < 4. Secondly, they are isotropic in their scal- 
ing. In contrast, in a bulk smectic the layer fluctuations 
within the layers scale differently than those along the 
normal to the layers. 

In H m —i the combination of relevant anharmonic 
terms and large disorder-induced u s fluctuations leads 
to strong anomalous elasticity By anomalous elas- 

ticity we mean that the full, anharmonic theory with con- 



stant K, B ss and A can, at small wavevector k <C 
(where £jvx is a non- universal length determined by ma- 
terial parameters, e.g., aerogel density), be effectively re- 
placed by a harmonic theory with wavevector dependent 
K, B ss and A. These are given by 

K(k) = Kk-^fK (k h /ki h ,k s /k^) , (3a) 
B ss (k) = B ss W z B f B (k h /k<?,k./k£) , (3b) 
A(k) = Ak-^f A (k h /ki\k s /ki>) . (3c) 

Bh s is not significantly renormalizcd; that is B^ s (k) = 
Bh Sl independent of wavevector. Here the anisotropy ex- 
ponents Cs = 2—(r/B + r/K)/2 and £/, = 2 — rj K /2. The ex- 
ponents, evaluated using the renormalization group and a 
high precision e-expansion were found to be: r\K = 0.50, 
n B = 0.26, and r/ A = 0.13. §,[§ We also predict that 
the anomalous exponents will obey the following exact 
scaling relation: 



1 + VA = W 2 + 2r )K- 



(4) 



The translational order of the system along the soft direc- 
tion is short-ranged and is characterized by the algebraic 
and anisotropic divergence of u s correlations, 



r z > r Slh 



<(<Mr)) 2 >= (lB(§(lB 1/2 ) Xa > r »> 



(#:)(§M#r) 1/2 ) Xh : r h >r. 



(5a) 



where we have defined Su s (r) = u s (r) 
1 - Vk + W 2 + Va = Vb + Vk, Xs,h = 



,1/2 



- u s (0), Xz = 
Xz/Cs,h , and 

£ 2 = K 2 B 1 J h 2 /(ABU 2 ). In the second equality for % z , 
we have used the exact scaling relation between the 77's. 
The exact scaling relation EqJJ, could be experimentally 
tested by using the more general expressions for X zi Xs: 
Xh in terms of all three exponents, and verifying that 
Vb , t\k , andi]A obey Eq.|]. Our e-expansion results for 
the 77's imply X z = 0.76 ,Xs = 0.47 and X h = 0.43. The 
fluctuations given in Eq.[5a|, like those along the hard di- 
rection, are disorder-, rather than thermally-driven. 

Despite this lack of translational order, our detailed 
calculations p] indicate that dislocation loops remain 
bound for weak disorder, and therefore the low tempera- 
ture phase replacing the columnar phase must be distinct 
from the smectic and hexatic, separated from them by 
a thermodynamically sharp dislocation unbinding phase 
transition. 

The stability of this exotic glass phase is contin- 
gent upon our assumption of long-ranged orientational 
order. We validate this assumption by calculating 
(|n(r) - n(0)| 2 ) = (|<9 z u(r) - d 2 u(0)| 2 ) and showing 
that it does not diverge as r — * 00 H . Although equili- 
bration into the ground state might be slow and therefore 
require field alignment, this orientational order would al- 
low experimentalists to investigate single domain samples 
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of HCBG. The anisotropic scaling information which is 
usually lost in a powder averaged x-ray scattering exper- 
iment would be retained in a single domain experiment 
allowing detailed tests of our predictions for t]k, t\b and 
VA- 

The scattering pattern in the _L plane, obtained from 
a single domain sample, would consist of a set of spots 
rather than the set of rings that one would expect from 
a powder sample. This pattern depends crucially on the 
relative orientations within the _L plane of the reciprocal 
lattice and the axis of stretch, e/j, which could vary from 
discotic to discotic since it depends on the microscopic 
interactions between the discs and strands. The intensity 
of a Bragg spot at a reciprocal lattice vector G 



J(G)« jTeKp[-<[G.(u(r)-u(0))P)/2] 



(6a) 



oc 



jf«*p{-[G£((u h (i 



u h (0)) 2 



G 2 s ((u s (r)-u s (0)) 2 )\/2} 



(6b) 

Unless G s = 0, the algebraically diverging u s fluctua- 
tions dominate the logarithmically diverging Uh fluctua- 
tions and the integrand is exponentially damped, leading 
to an anisotropically broadened Bragg peak. If however, 
G s = then the exponential becomes r -0 55 ™ , where 
n = Gh/Goh, with Goh being the magnitude of the small- 
est G lying on the hard axis, and for n < 3 the integral 
diverges as r — ► oo, leading to quasi-sharp peaks for those 
n's. We therefore predict two classes of hybrid columnar 
Bragg glasses. The first, which we call a commensurate 
HCBG, has some reciprocal lattice vectors that lie along 
the hard axis, and will exhibit a scattering pattern with 
peaks lying on the hard axis, with the first two quasi- 
sharp. In the second, incommensurate HCBG, class, all 
the peaks lie off the hard axis, and are anisotropically 
broadened by the contribution from the u s fluctuations. 
The smectic-like scattering pattern (Fig.^J), with a quasi- 
sharp peak on the first ring, will therefore only be ob- 
served for commensurate HCBG's. 

The dependence of the anisotropically broadened peak 
widths on the bare elastic constants is the same for both 
classes of HCBG. Setting the u — u correlation functions, 
given in Eq.^a], equal to G~ 2 and solving for (r z — £^) _1 , 
( r s = C^) _1 an( l ( r h = £/f ) _1 gives the width of the peak 
at G along each of the z, s and h directions: 
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- 1 {G 8 2 K/B SS )^ 
- 2 {K/B ss y' 2 {G s 



'K/B s 



(K/B sh )^ 2 (G s 'K/B 6 



,)x- h 



(7a) 
(7b) 
(7c) 



The temperature dependence of h could be used to 
determine the exponents t\k, t\b and tja since the bulk 
K(T), B SS (T) and B sh (T) in Eq.0 have temperature de- 
pendences that can be extracted from data on bulk ma- 
terials. A more direct way to observe the anomalous elas- 
ticity would be a direct measurement of the u — u correla- 
tion function I(q) oc (|u s ((5q)| 2 ), which can be obtained 



H for large q (i.e., q's with at least one component bigger 
than the corresponding inverse x-ray correlation length 
quoted above) by looking at an intermediate regime in 
the "tails" of the broad x-ray scattering peaks. In those 
tails, (i.e., for q = G + 5q with > \5q a \ > 



for at least one Cartesian direction a 
can be shown that 



J(q)« 



A(5q)q 2 



(B ss (5q)Sq 2 +K(Sq)Sqi 



(h,s,z) M), it 



(8) 



BhMl) 2 



Hence, the q dependence of B ss (q), K(q) and A(q) in 
Eq.|| could be tested directly by a fit of scattering data 
to these tails. 

A related experimental approach, which has the advan- 
tage of not being restricted to wavevectors larger than the 
inverse x-ray correlation lengths, but can, rather, explore 
arbitrarily small q's, is light scattering, which measures 
director fluctuations. These can be related to the u—u 
correlations via our condition <5n as d z u. This yields 



(l^(q)l 1 



(\Sn h (qW 



A(q)^ 



(B ss (cOq 2 +K( q )qi + B hsq 2 y ' 

(C(3)/2G 2 h q 2 Jq 3 , 
\k B Tq 2 G(q) + A,(q)qf G(q) 2 , 



(9a) 



commen. 
incommen. 
(9b) 



where A^(q) in Eq. 9b is renormalized, q depen- 
dent tilt disorder variance obeying the scaling law 
A„(q) = A h q-^ + ^f Ah (q h /q<»,q s /qi>) and G(q) - 
l/(B s hq 2 + "fq 2 + Bhh<lh)- The commensurate and in- 
commensurate cases differ because in the commensurate 
case there is a random field acting on Uh, while in the in- 
commensurate case there is no random field, leaving the 
random tilt as the dominant disorder. 
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